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J> . Abstract 

' The interconnection between quantum mechanics and probabilistic classical me- 

chanics for a free relativistic particle is derived in terms of Wigner functions (WF) 
for both Dirac and Klein-Gordon (K-G) equations. Construction of WF is achieved 
by first defining a bilocal 4-current and then taking its Fourier transform w.r.t. the 
relative 4-coordinate. The K-G and Proca cases also lend themselves to a closely 
parallel treatment provided the Kemmer-Duffin /3-matrix formalism is employed for 
the former. Calculation of WF is carried out in a Lorentz-covariant fashion by stan- 
dard 'trace' techniques. The results are compared with a recent derivation due to 
J^Z Bosanac. 
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1 Introduction: Wigner Function As Link Between 
QM and CM 

Since the advent of Quantum Mechanics (QM), there have been several investigations 
aimed at understanding Schroedinger equations and their relativistic counterparts from 
stochastic /statistical premises in classical mechanics, epitomised by Focker-Planck [1,2] 
and Liouville equations [3-6]. In the process most such investigations have had to resort at 
some stage to some sort of ansatz designed to introduce the 'square-root' of a probability 
function [1-4], or make a more direct use of the wave function from the outset [5,6]. This 
of course does not address the core issue of the foundations of quantum physics, viz., the 
existence of complex probability amplitudes obeying the "superposition principle", and 
the need to invoke the "collapse hypothesis" to explain the outcome of a quantum mea- 
surement. Even without any commitment to resolve this central mystery of conventional 
quantum physics, it is both interesting and illuminating to explore the classical roots of 
quantum mechanics by comparing and relating their respective mathematical structures 
as in several contemporary investigations [3-6]. 

Among the several approaches to an alternative 'derivation' of QM, we would like 
to pick out for explicit scrutiny, a recent formulation due to Bosanac [4] in which the 
author makes essential use of the probabilistic (instead of 'trajectory-oriented') approach 
to Classical Mechanics (CM) via the Liouville equation, wherein the Uncertainty Principle 
is sought to be incorporated at the outset, as a means of capturing a good chunk of QM. 
Our purpose in this regard is more in the nature of an 'inverse' (to [4]) approach, which 
stems from a desire to give a more direct role to the QM wave function born out of its 
'being already there', thereby according it a more open role in the construction of the 
Wigner Function [7], as in other orthodox approaches [5,6]. As we shall see, this approach 
has the advantage of presenting the treatment of different types of relativistic equations 
(Dirac, Klein-Gordon, Proca, etc) within a unified framework. 

The central quantity that satisfies the Liouville equation is the Wigner-function [7] 
which in itself was conceived within the QM framework, as a bilinear in the 'wave func- 
tion'. Granting the existence of the wave function in turn implies the existence of its 
Fourier transform which represents its momentum space counterpart, and the two to- 
gether guarantee a built-in effect of the Uncertainty Principle, as clearly recognized in 
[4]. Using such a (hybrid) approach in which the prior existence of the wave function 
is implicit in the structure of the Wigner function [7], Bosanac [4] claims to derive the 
Schroedinger equation from the Liouville equation, consistently with current conserva- 
tion, apparently as an inverse(?) problem. An extension to the corresponding relativistic 
problem, rather surprisingly, seemed to yield only the Dirac equation (spinor), but not 
the Klein-Gordon (K-G) equation (scalar), a result whose mystery we seek to clarify in 
this paper and (hopefully) pave the way to a generalization (see below). 

Specifically, we have been actuated by a desire to re-examine the issue of Dirac versus 
K-G [4], but with a more conscious recognition of an active role of the 'wave function' 
at the outset, and in the process we find that the Liouville and current conservation 
equations for both the Dirac and K-G cases follow from identical premises, provided one 
has no qualms about the free use of multi- component wave functions in both Dirac and 
K-G cases, and not (selectively) for the Dirac case only [4]. Indeed it is easily seen 
from a detailed perusal of the steps in [4] that the author's success in deriving the Dirac 
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equation is crucial to his ansatz of a 4-component spinor, without which the necessary- 
consistency checks would not be forth-coming. On the other hand he did not seem to 
consider a similar possibility for the K-G case, viz., the existence of the less familiar 
Kemmer-Duffin equations for integral spins [8], wherein the multicomponent structure 
of the corresponding wave function satisfying a first order equation as in the Dirac case 
seems to be the right form for putting the K-G and Proca cases on parallel footing with 
Dirac. Of course the question of whether, inspite of a multicomponent ansatz for the K-G 
wave function, a bosanac-like treatment [4] of working backwards from the Liouville-cum- 
current-conservation equations to the K-G equation would have succeeded, is still open, 
since the Kemmer-Duffin /^-matrices are singular, unlike the Dirac matrices which are not 
(see Sec. 4). 

In carrying out this exercise we have been guided solely by the orthodox (Dyson [9]) 
philosophy that a relativistic quantization of a single particle has inherent problems of 
consistency, whose resolution hinges on the simultaneous existence of an infinite number 
of particles, i.e., on a field concept. Indeed, Dirac's postulation of the infinite sea (of 
negative energy particles) in the vacuum was enough ground for Pauli-Weisskopf [10] to 
reinterpret the single particle density from the K-G equation as an average charge density 
(which could 'locally' have either sign) [9]. Of course from a practical point of view, a 
single particle interpretation has considerable appeal, but it requires the 'protective cover' 
of Feynman's positron theory [11] which mandates ± time directions to be associated with 
± energy propagations respectively. 

In what follows, we shall keep in the forefront, the results of [4] for pointwise compari- 
son with the derivation to be presented here under the orthodox 'field' paradigm [9] which 
is routinely available in any standard text [12]. In Sec. 2, we give a 'direct' derivation of 
the Wigner function for the Dirac equation and thence the corresponding Liouville plus 
current conservation equations, a procedure whose steps are roughly in 'inverse' order to 
those of [4]. In Sec. 3, we do the corresponding exercise for the K-G and Proca cases via 
the Kemmer-Duffin [8] formulation [13], and in so doing, bring out the strong similarity 
of the two systems. Sec. 4 concludes with a short discussion, including a comparison with 
[4]- 



2 Wigner Function for the Dirac Equation 

We begin by writing the Dirac equation for a free field ip(x) in the 'Euclidean' covariant 
notation as [13] 

(7A + ™)#£) = 0; i;(x)(m-d^) = 0. (2.1) 
Now define a bilocal current j^,(xi,x 2 ) as 

jn(xi,x 2 ) = i>(x 1 )i-/^(x 2 ) (2.2) 
Next put Xi = x — q and x 2 = x + q as in [4], and define a 4- vector V^(x,p) as 

V^(x,p) = ir- A J d A qe 2iM 4>(x + q)i^ip(x - q). (2.3) 

The fourth component of may be identified with the standard WF p. As a check on 
this quantity, its integral w.r.t. d A p yields the usual 4-current as 

J^x) = / d 4 pV^(x,p) = i>(x)i-f^(x) (2.4) 
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from which the standard current conservation equation is recovered as 



dMx) = 0. (2.5) 

Next we make a Gordon reduction of (always possible for a free field i>(x)) so as to 
bring out the separate contributions and V* arising from the 'connective' and 'spin' 
currents respectively: 

V'{x,p) = 7T~ 4 f d 4 qe 2ip -y(x + q)^{x - q); (2.6) 

V'(x,p) = J d 4 qe 2i ^(x + q)^(x - q). (2.7) 

We check from these that two separate conservation equations follow: 

d lt Vfa,p) = d li VZ(x,p)=0, (2.8) 

the former being a consequence of the K-G equation for each component of ip(x) and 
■0(x), and the latter being identically satisfied. 

We now seek to use the Gordon reduction as a key element for deriving the Liouville 
equation a la [4], and in so doing bring out its role in putting our procedure in closer 
correspondence to [4]. For, the central issue is one of demonstrating that the density 
function p(x,p) and the current function Vi(x,p) are related by the 'classical' velocity 
factor vi, so as to ensure the 'correct' structure of the Liouville equation [4]. To that end 
we note from (2.3,6,7) that the effect of Gordon reduction in (3) may be expressed in 
operator form sandwiched between ip and ip: 

il» = VjL + ^-0< X) <V (2.9) 
m 2m 

The spatial part of this equation is seen as the precise counterpart of eq.(27) of [4] (after 
integration w.r.t. d 4 pd A q), with the following identifications: The convective term p^ in 
(2.9) has a space part which corresponds to the j% term of [4], since it bears the ratio 
to the time component po, as the correct factor connecting the current function Vi(x,p) 
to the density function Vo(x,p) . The space part of the full 7 M term collectively represents 
the two terms aid and a^e in eq.(27) of [4]. Finally the last (derivative) term in (2.9) 
corresponds to the V x $ term in eq.(27) of [4] since it has a zero divergence by itself, 
as seen from the second part of eq.(2.8), and therefore has no separate influence on the 
equation of continuity. Therefore, following the logic of [4], this term has no effect on the 
anatomy of the Liouville equation which is thus entirely contained in the first (convective) 
part of eq.(2.8). To reveal this structure more explicitly, the p(x,p) function of [4] reads 
in our 4x4 Dirac matrix notation, as follows: 

p(x,p) = V (x,p) = J d 4 qe 2ip - q ip(x + q)^{x - q). (2.10) 

which satisfies the Liouville equation 

d p(x,p) + Vidip(x,p) = (2.11) 

that is merely a paraphrase of the first part of eq.(2.8) representing the conservation of 
the convective part of the function. 
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To end this section, we outline a construction for the p function analogously to eqs.(51- 
53) of [4], using a fully covariant procedure. To that end, we write down Fourier repre- 
sentations for the Dirac spinors on the lines of [12], but adapted to the Euclidean phase 
convention [13]: 

^(xi) = (2tt)- 3 / 2 | ^^/^il^^n^lPiJe 111 " 1 + rf; m " ri (pi)e- rai '1; (2.12) 

$(x 2 ) = (2tt)- 3 / 2 | d 3 p 2 (m/E 2 y/% 2r2 u r 2(p 2 )e- ix ^ + d P2r2 v r *(p 2 ) x e™]; (2.13) 

where b and d are (c-number) electron and positron amplitudes respectively and u and v 
are the corresponding positive and negative energy spinors. The latter are more conve- 
niently expressed in terms of Lorentz boost operators as follows [12]: 



u^( Pl ) = [(m - i 7 .p 1 )/(2 v /m(m + E l ))}u r ^ (0); (2.14) 



v ri (Pi) = [(m + i>y.p 1 )/(2yjm(m + E^y^O). (2.15) 

These Lorentz boost factors are proportional to energy projection operators for the cor- 
responding spinors. To specify their spin states, we can further multiply these quantities 
by the spin projection operators P n [12] 

P n = [l + i 7 .n 75 ]/2, (2.16) 

where n is a unit 4- vector representing the direction of spin [12]. After these manipulations 
on the spinors through the sequence of equations (2.10-16), and their substitution in the 
defining eq.(2.3) for V^, the spinor dependence can be totally removed by 'tracing'. To 
simplify these expressions note that only the b*b and d*d terms will survive, while the 
cross terms will drop out of the traces, since the ± energy spinors at a given momentum 
(zero) are orthogonal to each other. The rest of the simplification is routine and the result 
can be compactly expressed in terms of two quantities Tr^ defined as follows: 

Tr ± = ±m(pi + p 2 ) M + ^ApoPuPgX ^ 
" 2yjE 1 (m + E 1 )E 2 (m + E 2 ) 

Before substituting these quantities in eq.(2.3) for V^(x,p), note that of the two terms 
on the right, the first (convective) term has the desired structure for the 'Liouville form', 
eq.(2.11), but the second (spin) term has not. However, in accordance with the discussion 
immediately following eq.(2.9), in line with [4], one need not insist on the 'Liouville 
connection' separately for both the terms since the 'spin' term (proportional to n) can be 
traced back to the 'derivative' term of (2.9) which has no separate effect on the equation 
of continuity. Therefore we need only the p — 4 component in eq.(2.17) above to extract 
the p-function fully, whence the corresponding 'current' function V$ can be generated by 
multiplying with Vi = Pi/po only. 

To proceed further, integration over d 3 pi and d 3 p 2 may be simplified via the transfor- 
mations pi t 2 — p' ± k/2, remembering that x 2>1 — x ± q. Further, the 3D integrals may 
be converted to 4D integrations over p' and k via the identity 

d 3 Pl d 3 p 2 = 2E 1 2E 2 d A p'd' l k5{p' 2 + m 2 + k 2 / A)5{2p' .k) (2.18) 
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making use of the result 5 (a) 5(b) = 5((a + b)/2)5(a — b). Integration over d A q now gives 
& A (p' T p) for the terms b*b and d*d respectively; subsequent integration over d 4 p' gives 
p' = ±p. Collecting all these results and substituting in (2.10), gives 

p(x, p) = 7T- 3 J d A k5(p 2 + m 2 + k 2 /4)5(2p.k) 

x[b*(p + k/2)6(p - k/2)e lk - x + d*(p + k/2)d(p - k/2)e~ lk - x ] 

x(2 m p + 2in.(p x ^ {m + ^Z + E 2 ) < 2 " 19 > 

This expression corresponds to eqs. (51-53) of [4], after taking due account of the difference 
in the choice of the respective basis functions (Dirac- vs. Pauli- types), as well as the 
covariant vs. non-covariant notations. Note that the integration in (2.19) is basically 
angular dQk in content, like in [4], since the 4D integration d A k is also multiplied by the 
two ^-functions which take care of the magnitudes of both k and k . Our 3- vector k is 
2k of [4], and n is u> of [4]. Finally, our expression for p compactly covers both ± energy 
components. 

We next turn to the K-G and Proca cases via Kemmer-Duffin equations. 



3 Wigner Function for Integral Spin 

Because of the relative lack of exposure of the Kemmer-Duffin equations [8] in the contem- 
porary literature, it will first be useful to recall some basic results which are summarized 
in Roman [13] in the Euclidean phase convention, from which we shall freely draw on 
several results without much comment. The first order differential equation for integral 
spin is expressible in (5- matrix notation as [13] 

(P II d II + m)^(x) = 0; $(m-dM = (3.1) 

where ip(x) is a multicomponent field, each of whose components ipi satisfies the K-G 
equation, and ip = ■0^4 and ^4 = 2/3| — 1. The spin-0 and spin-1 cases correspond to 5- 
and 10- component irreducible representations respectively. These components are 

Spin — : ipi(i — 1 — 4) = — <9 M 0(x)/m; ip 5 = 

5^71-1:^(1 = 1-6) = -F^/m- ^(i = 1 - 4) = M (3.2) 

The various quantities like bilocal 4-current j /J ,(x 1 , x 2 ), phase space 4-current V^(x,p), and 
local 4-current J^(x) follow very closely the corresponding constructions (1.2-5) for the 
Dirac case, with the replacements: 7 M by (5^ ; 74 by 774. Thus the phase space 4-current is 
given by 

V„(x, p) = vr" 4 J d i qe 2ip - q ^(x + q)i(3^(x - q) (3.3) 

and its break-up into convective and spin components is achieved via 'Gordon reduction' 
analogously to the Dirac case. However the algebra of the f3 matrices is a bit more involved 
than that of the 7 matrices because the former are singular. Thus [8,13] 

(3xMv + (3„Mx = (3x5^ + (3„5»x; r,l = l (3.4) 
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The 'spin' operator S^ v is given by 



iSp, = PvPu - PuPv (3.5) 

and, except for a factor of half, has exactly the same role as that of a^, including 
its antisymmetry property. For writing down the counterpart of Gordon reduction, the 
following equations are needed: 

(fy - PuPM^x) = 0; <Kz)(c^ - frPM = (3.6) 

With the help of (3.1) and (3.6), and making use of the definition (3.5) for the 'spin', the 
matrix (5^ in (3.3) may be written in a form analogously to (2.9) as 

t(3, = ^ + d^^ (3.7) 
m m 

which shows a clear division between the convective and spin parts, with the latter having a 
'derivative' structure whose divergence vanishes by itself, as in the second part of eq.(2.8). 
Therefore this term has no direct role in the Liouville equation like in the Dirac case, 
exactly as in [4]. 

So far this /3-formalism for integral spin covers both spin-0 (K-G) and spin-1 (Proca) 
cases, and one may wonder why for the former there should be a spin term S^ v at all, 
in the Gordon reduction (3.7). Actually this term is effectively zero for the spin-0 case 
since, as has been explained in Roman [13], for the 5x5 representation the eigenvalues 
of S3 = S12 are non-measurable. Indeed the probability density [13] 

p = -itytfc + = -i - <f>*d (f>)/m (3.8) 

vanishes in the eigenstates that correspond to S3 = ±1, so that the presence of the 
S/jus-term in the spin-0 case is purely symbolic. Substituting this information in (3.3), 
and noting from (3.8) that only the components ^4,5 survive for the phase space density 
function p(x,p) for the K-G case, the final result works out as 

p(x,p) = —n~H J d A qe 2ipq [d Q (t>*{x + q)<j>(x - q) - (f)*(x + q)d (f>(x - q)} (3.9) 

while eq.(3.7) ensures that the phase space current function V(x,p) is related to p(x,p) 
by the factor p/po = v, so that the structure of the Liouville equation is preserved. As a 
further check, the integration of (3.9) over d 4 p reproduces the K-G density function p(x) 
of (3.8) exactly. The last quantity is of course not positive definite, as expected [13] for 
the K-G case. Again as in the Dirac case, it is possible to write down an expression for 
p(x,p) closely analogous to (2.19), except for the absence of the 'spin' term, in the third 
line. We omit this routine step for brevity. 

For the Proca case, the 'spin' term in the Gordon reduction (3.7) will now contribute to 
the p(x, p) function, while the Liouville connection is again maintained via the convective 
Pn term only. From the structure of p(x,p), which corresponds to p — 4, it is seen 
quite clearly that the 'spin' contribution arising from the second term of eq.(3.7) in the 
Proca case is proportional to S.E x H, where E and H are the 'electric' and 'magnetic' 
components of F^ v respectively. Again we skip the Fourier representation (2.19) for the 
full p(x, p) function, as this step will not illuminate the structure any further beyond the 
above demonstration of a close parallelism between the spin- 1/2 and integral spin (0,1) 
cases. 
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4 Discussion and Conclusion 



This study was motivated by a desire to understand some results due to Bosanac [4] 
in which he claims to derive the Dirac equation, starting from the (classical) premises 
of the Liouville equation. What intrigued us was the assertion [4] that the 'derivation' 
leads to the Dirac equation, instead of to the Klein-Gordon equation. To examine this 
question more closely, we were led to start from the opposite direction to [4], viz., a 
'direct' derivation of the Liouville equation through an explicit construction of the Wigner 
Functions [7] themselves, as in more orthodox approaches [5,6], with a view to throwing 
some light on this apparent asymmetry between the Dirac and the K-G cases. Rather 
surprisingly, our result indicates a close parallelism between both spin-1/2 and spin-(0,l) 
cases, brought about by the Kemmer-Duffin [8] formalism [13] for integral spin as a natural 
counterpart to the Dirac equation for spin-1/2. In both cases, we have made essential 
use of the Gordon reduction into 'connective' and 'spin' contributions to the 4-current, 
wherein the convective part gives the desired 'Liouville ratio' po : p between the density 
and current functions in phase space, but the spin part does not. However the latter 
has a derivative structure by virtue of which its 4-divergence vanishes identically, so as to 
yield a 'spin' current conservation by itself [eq.(2.8)], and hence has no direct effect on the 
dynamics of the Liouville equation as such, which is entirely governed by the convective 
term only. This result seems to accord rather well with the corresponding result of Bosanac 
[4], viz., his eq.(27), despite a vast difference in the respective procedures. In particular, 
since Gordon reduction in our derivation has proved crucial for a neat identification of the 
origin of these two terms as convective and spin contributions respectively, it may be safely 
assumed that the same mechanism underlies the derivation in [4] as well, irrespective of 
the manner of the derivation. 

Without going into the philosophical aspects of the classical foundations of Quantum 
Mechanics, it is pertinent to ask why only the Dirac equation has been amenable to 
the Bosanac derivation [4]. For, it will be rather naive to argue that no substantial 
quantum mechanical input has gone into the latter approach, since the very structure of 
the Wigner Function as a bilinear in the (quantum mechanical) wave function is a tacit 
acknowledgement of the existence of the latter, with the Uncertainty Principle as a built- 
in corollary of the same. Further, the structure of the Dirac equation is implicit in the 
ansatz [4], eq.(16), in terms of Pauli spinors Fi, F2, with anticommutation relations (23-25) 
characteristic of Dirac matrices in the Pauli 2-component form, all of which ingredients 
have been crucial to the (indirect) derivation of the main eq.(27). Indeed the entire 
exercise [4] is a fine example of how to work backwards from the (bilinear) Liouville 
equation to the (linear) Dirac equation, albeit in the Pauli 2-component form. Yet it goes 
without saying that the ansatz of a 4-component spinor (satisfying a Dirac-type equation 
(16) [4]) has been very much of an input in the derivation. 

A comparison with our 'direct' derivation would suggest that a possible reason why 
the K-G equation was not amenable to the Bosanac derivation was the missing element 
of a multicomponent wave function in his formulation of the K-G problem, something we 
have sought to rectify through the Kemmer- Duffin /3-matrix formalism [8] which restores 
the formal analogy to the Dirac case. However, and this may be important, even if the 
Bosanac formalism had included a multicomponent wave function for integral (including 
zero) spin field, it is doubtful if working backwards from the Liouville to the K-G equation 
would have proved successful in this case, since the (5- matrices are singular, unlike the 
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Dirac matrices, so that a routine inversion of steps on the lines of [4] might well not have 
been possible; (we did not succeed in this effort). But there seems to be no problem in 
the 'forward' steps (from QM equation to Liouville) going through in the K-G and Proca 
cases with equal ease, as the above derivation indicates. 

The foregoing comparison would suggest that the structure of Quantum Mechanics is 
probably richer and more varied than can be fully captured in one go, through alternative 
formulations. And yet the basically indeterministic nature of QM poses a continuous 
challenge to innovative ideas designed to capture its full flavour from 'objective-realist' 
premises, so as to provide a more 'classical' approach to its uncanny predictions. Ref. [4] 
may be regarded as one more effort in this direction, as part of an ongoing movement 
[14] for resolving the foundational paradoxes of orthodox QM which has acquired fresh 
momentum with heavyweights for and against the Copenhagen Interpretation arrayed in 
comparable strengths on both sides of the Great Divide [15]. 
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